arXiv:1504.02733vl [hep-ph] 10 Apr 2015 


Prepared for submission to JHEP 


Lepton-pair production in ultraperipheral collisions at 
AFTER@LHC 


J.P. Lansberg“ L. Szymanowski 6 J. Wagner 6 

a IPNO, Universite Paris-Sud, CNRS/IN2P3, F-91f06, Orsay, France 
b National Centre for Nuclear Research (NCBJ), Hoza 69, 00-681, Warsaw, Poland 

E-mail: Jean-Philippe. Lansberg@in2p3. fr, Lech. SzymanowskiOfuw. edu.pi, 

Jakub.WagnerSncbj.gov.pi 

Abstract: We investigate the potentialities offered by the study of lepton-pair production 
in ultraperipheral collisions at a fixed-target experiment using the proton and ion LHC 
beams. In these collisions, exclusive or semi-exclusive lepton-pair production can be used 
as luminosity monitor as well as a check of the equivalent-photon approximation, via the 
measurement of the Bethe-Heitler cross section. It can also serve as a probe of the inner 
hadron structure via the measurement of the lepton-pair azimuthal asymmetry which is 
sensitive to the timelike virtual Compton scattering. We also briefly discuss the possibility 
offered by the study of ry c production. Finally, we outline the possibilities for lepton-pair 
production by Pomeron-Odderon fusion in exclusive pp and pA collisions. 
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1 Introduction. 

With the advent of RHIC and the LHC, the first experimental studies of Ultra-Peripheral 
Collisions (UPC) have successfully been carried out. The STAR collaboration first mea¬ 
sured the po production cross section in AuAu UPC [1], then they measured the di-electron 
production cross section in UPC [2], A little later the PHENIX collaboration released their 
first study of the <//?/; cross section as well as high-mass di-electron production [3]. More re¬ 
cently, further studies have been carried out at the LHC by the ALICE collaboration [4, 5] . 

Attempts to isolate UPC in Inin collisions in the fixed-target mode at the SPS have 
been made 1 . They were not conclusive, most probably because of the limited nucleon- 
nucleon centre-of-nrass (cnrs) energy, on the order of 20 GeV, resulting in typical photon- 
nucleon cnrs energies below 3 GeV. 

In this context, we investigate the possibility to study lepton-pair production in ultra¬ 
peripheral collisions at a fixed-target experiment using the proton and ion LHC beams [6] - 
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referred thereafter to as AFTER@LHC. In such collisions, one can investigate specific reac¬ 
tions where one of the colliding particle serves as a (coherent) emitter of a photon and the 
other serves as a target. Such photon-hadron collisions can be (semi-)exclusive, resulting, 
for instance, in the sole production of a lepton pair. This process can be used to monitor 
the experiment luminosity since it mainly comes from the Bethe-Heitler (BH) process [7], 
whose cross section is well-known. By looking at the target rapidity region, it can also be 
used to scan the domain of validity of the Equivalent-Photon Approximation (EPA). It can 
also be used to probe the inner structure of the target, through the interference between 
the BH process and the Timelike Compton Scattering (TCS) [8-10]. Such an interference 
- measurable via the analysis of the azimuthal anisotropy - indeed involves contributions 
from the Generalised Parton Distributions (GPD) [11-18]. 

The structure of this article is as follows. In section 2, we present the main charac¬ 
teristics of the UPCs and the corresponding photon fluxes in a fixed-target mode on the 
LHC beams. In section 3, we briefly discuss the cross sections for production of lepton 
pair via the BH process. In section 4, we discuss how the contribution from TCS can be 
extracted and how they can help to unravel information about the inner proton structure. 
In section 5, we briefly discuss the potential competing hadronic process resulting from 
photon-odderon fusion. Finally, we present our outlooks and conclusions. 

2 Ultraperipheral collisions in a high-energy fixed-target experiment 

2.1 Generalities on photon-induced reactions in ultraperipheral collisions 

Relativistically moving charged hadrons are accompanied by electromagnetic fields which 
can effectively be used as quasi-real-photon beams. At very high energies, these photons 
are energetic enough to initiate hard reactions, just as in lepton-proton colliders. 

The virtuality, q 2 = —Q 2 , of these photons is small, Q 2 < 1 /R 2 , where R is the radius 
of the charge. More precisely, Q < 0.28 GeV for protons (R m 0.7 fm) and Q < 0.06 GeV for 
nuclei (Ra ~ 1.2 A 1 / 3 fm) with a mass number A > 16. These photons, which are emitted 
coherently, are almost on mass shell, and their emission can be theoretically treated in the 
EPA (see e.g. [19]). 

Seen from a target at rest, the energy of these photons can become significant if the 
energy of the moving charge, i.e. the beam energy, becomes ultra-relativistic, as at the 
LHC. At rest, the coherent photon cloud of an heavy ion is on the order of 30 MeV. 
Boosted a few thousand times (qpb — 2940), these photons have an energy close to 100 
GeV in the laboratory frame. It is of course much less than what can be achieved at the 
LHC in the collider mode, but close to the experimental condition at RHIC with colliding 
beams of 100 GeV. It is anyhow enough to produce hard dileptons as well as vector mesons. 

The energy spectrum of these photons depends on the boost with respect to the ob¬ 
server as well as the impact parameter b - it is understood that the observer or the probe 
is outside the charge distribution. Using the EPA method, one gets [20] that the flux as 
function of the photon momentum fc, of b and 7 (the Lorentz factor of the hadron - or 
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nucleus - in the frame where k is measured) reads 
dn Z 2 a era uj(b,k) 2 r 


dkd 2 b 


7 T 2 kb 2 


Kf(u(b,k)) + ^K 2 (u;(b,k)) 


( 2 . 1 ) 


where a e m is the QED coupling, Z is the nucleus charge, oj( b,k) = kb /7 and K\p are 
modified Bessel functions of the second kind. 


Table 1: Relevant parameters for AB UPCs at AFTER@LHC, at RHIC and at SPS: (i) 
nucleon-nucleon cms, t/siviv (h) luminosity, Cab, (iii-iv) colliding hadron energies, E l j/ h B , 
in the laboratory frame, (v) Lorentz factor between the colliding hadron rest frame and 
cms, 7 = ^snn /(2 mjv), (vi) Lorentz factor between both colliding hadron rest frames, 
7 = s NN /(2m 2 N ), (vii-viii) inverse of the colliding hadron effective radii (giving the typical 
photon cloud energy in the emitter rest frame) (ix-x) photon “cutoff energy” in the target 
(resp. projectile) rest frame, ElJ (resp. E/// r f^) (xi-xii) “maximum” photon-nucleon 
cms energy where A (resp. B ) is the photon emitter, (resp. (xiii) photon 

“cutoff energy” in the cms, with both A and B emitting a photon coherently (xiv) 

“maximum” photon-photon cms, 
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100 

115 
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16 
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15 

30 
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2.8 x 10~ 3 
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25 

1.5 
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Inin 

n/av 

17 

n/av 

160 

m N 

9.23 

170 

16.9 

2.87 

2.4 

0.16 

0.31 

PbPb 

n/av 

17 

n/av 

160 

m N 

9.23 

170 

13.9 

2.36 

2.1 

0.13 

0.26 


Unless b is smaller 3 than R , the strong suppression of the flux by the Bessel functions 
when oj gets of the order of unity implies that k should be smaller than inverse radius times 
the Lorentz factor 7 . The larger 7 is and the smaller the emitter is, the harder the flux 
is. We also note that the energy spectrum kdn/dk is constant for kj 7 fixed. To fix the 


idea, one usually considers a maximum photon momentum /t max — 


he 


below which the 


Remitter 

emission are likely coherent and therefore characterised by a flux proportional to Z 2 . This 


2 For Arp and Op luminosity with AFTER@LHC, we conservatively assumed the same extracted flux of 

Ar and O as for Pb, i.e. 2 x 10 5 Pb/s. See also [6, 21-23] 

^Otherwise, (i) one cannot consider the entire nucleus charge Z, (ii) hadronic interactions may be more 

important that the photon-induced ones, (iii) the probably for the colliding objects to break-up may also 
be important. 
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quantity should not thus be considered as a sharp cut-off above which photon emissions 
are forbidden. If the photon is considered in the hadron cms, one has 7 = y/sNN /(2mjv)(= 
7 cms ). If the photon is considered in the target B rest frame, 7 = s n n / {2m 2 N )(= r y A ^ B ). 
Table 1 summarises the relevant parameters characterising ultra-peripheral collisions at 
AFTER@LHC, at RHIC and a SPS in fixed-target mode. 

For comparisons with photon-induced reactions in the more conventional lepton-hadron 
collisions, it is usually more instructive to look at the maximum of the 7 N cms en¬ 
ergy) y/s™jY X - To do so, we “boost ” 4 fc max in the “target” nucleon rest frame, to obtain 

Ey nfax = 7 A ^ B k max , where y/s^ = y/2E^ a ^ st mN = \/sNNk max /mN. We note that the 
photon-energy “cut-off” obtained with these dimensional arguments for Pbp, i.e. 74 GeV, 
is remarkably close to the peak in the energy spectrum obtained using a more realistic 
model in a recent study of the Bremsstrahlung spectrum of ions in AFTER@LHC [24], i.e. 
80 GeV. 


2.2 Photon fluxes 

Taking into account the smallest possible impact parameter for a given colliding system, 
pp. pA or AB, as well as the charge distribution through a form factor in the proton case, 
one obtains different formulae for the flux integrated in b. In fact, the photon fluxes do not 
formally factorise since b m i n depends on the radius of both colliding objects, except in ep 
collisions, where one can reasonably neglect the electron radius. 

Along these lines, one should normally have for pp collisions, 6 m i n ~ 2 x R p ; for pA 
collisions, b m ; n ~ R p + Ra ; and for AB collisions b m i n — Ra + Rb- Whereas it is acceptable 
to approximate R p + Ra to Ra, it does not seem justifiable to use Rpb for PbPb collisions, 
for instance. In addition, in pA collisions, it is also problematic to use a different 6 m ; n 
when one considers the proton emission or the ion emission. In both cases, one should use 
R p + Ra, or perhaps R,\. 

Integrating Eq. (2.1) over b, one has [20] 


dn 

dk 


2Z 2 a em 

irk 


k’(frmin) k~)K q (cu( 6 m i n , k)^I\\ [ui(b m i n , k)) 


w(^min) k) 


(Al (cu(6 m i n , ^)) Rq k)') 


( 2 . 2 ) 


To avoid any confusion with the choice of the frame, it is useful to work with the 
momentum fraction or light cone coordinate, = k/ph — k/[pfM^ r), where ph is the 
momentum of the hadron emitting the photon. One trivially obtains 


dn 

dx 1 


k dn 
rc 7 dir —.7)". 


(2.3) 


4 In fact, the procedure looks more as if the emitter is Lorentz contracted as I/ 7 , rather than the photon 
momentum boosted. The results are however similar once one considers the emitted photon as slightly 
off-shell, with a momentum, in the rest frame of the emitter, as (fc max , 0 , 0 , 0 ). 
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The relation between the (differential) hadron-hadron cross section, ( d)a hAflB , and the 
(differential) cross section for a photo-hadron scattering (h^ or hs), (d)a' yhA ’ B , is naturally 
given by the following convolution with the photon flux 


da hAhB = J d^(^da^( S7hfl (fc y )) + ^d^(5 7fcA (fe y )) 


dn hB 


v 7 

= Jdx 1 (^da^ B M k 1 ^)) + d ^- da' )hA {s lhA {k 1 {x 1 ))) S j. 
By analogy with the parton model formulae, one can thus write: 

^ = 


(2.4) 


(2.5) 


and interpret the latter as an equivalent-photon PDFs off the hadron (or the ions) by 
leaving the emitter intact. 



x i 


Figure 1: ^ for photons from protons (blue dot for Drees-Zeppenfeld, green dash for 
bmin = Rp and red solid for b rn \ n = 2 R p ) and from Pb (orange double dot dash for b m in = 
-Rpb, light blue dot dash for 6 min = Rpb + Rp and black double dot for 6 m ; n = 2i?pb ) 
divided by Zp h . These are compared to the gluon (thin pink long-dash, divided by 100) 
and u quark (thin purple dot-dash) MSTW PDF in the proton. 


In the case where the emitter is a proton, one could think that it is more accurate 
to take into account the spatial distribution of the charges through a form factor, which 
leads [25] to 


dn 

dxry 


a 


2nx~ 


(l + (1 — x 7 ) 2 ) 


11 3 3 1 

111 C — - + — — -77 -77 

6 c 2 c 2 3c 2 


( 2 . 6 ) 
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where c = 1 + °'q 2 GcV , and Qg = M 2 x 2 . This is only strictly relevant for ep collisions and 
at rather large x 7 . In pp or pA collisions, one indeed needs to take into account the radius 
of the other colliding object. The impact parameter b is therefore not getting close to R p . 
We will refer to this choice by “Drees-Zeppenfeld”. 

Fig. 1 shows a comparison (up to the charge factor Z 2 ) between the equivalent flux 
from a proton and from a lead ion. For < 10~ 3 , the spectra show a similar behaviour, 
with a magnitude differing by less than a factor 2 . At larger x 7 , the proton spectrum is 
clearly harder because of the smaller proton size. This is the expected behaviour. 

One could also consider the case where the proton emitting the photon breaks apart. 
In this situation, the photon is in fact radiated by the quarks and the corresponding photon 
PDF can then approximated by [26] 


hreak-up^ q2) = 


Q 


/ — p 7 q ( x /y) q(y, Q 2 ) + Q 2 ) 

271 QbaJxy 1 


dy 


(2.7) 


with P iq (z) = e 2 (1 + (1 — z) 2 ) / z, Qq is an energy cut-off, and q(x, Q 2 ) is the quark PDFs 
in the proton. We will not consider this possibility further in this study, although such a 
process could contribute to the semi-exclusive case with still a large rapidity gap between 
the dilepton and the proton emitter. Another possibility, which we will discuss in section 5 
is to have an elastic hadronic reaction via the exchange of a pomeron (or possibly an 
odderon). 


2.3 Fluxes and the rapidity dependence of the produced particles 
2.3.1 Single-photon case 

If we consider the total cross section, cr lh , to photo-produce a particle of mass Q, we note 
that it can only be function of s^h and Q 2 since it is already integrated on the final state 
variables. However, nothing prevents us to perform a change of variable in Eq. (2.4) from 
fc 7 to a final state variable of the 7 -hadron process, for instance the rapidity y of the 
produced particle, keeping the other fixed 5 6 . By momentum conservation, y would enter 5 
via fc 7 (y) in s^h 

Indeed, in our case, it is clearly instructive, in order to understand where the produced 
particle by photon-induced collisions fly, to consider the flux as a function of the particle 
rapidity y in the laboratory frame. For instance, we anticipate [ 6 ] that rapidities (in the 
laboratory frame) from -4 to +1 should be accessible. 

For fixed qx and Q 2 , /c 7 (y) is also fixed. From Eq. (2.4), one trivially obtains 

dfc 7 \ da lhB {s^ hB {k^)) + dcr 7 ftA (s 7 hA (A: 7 )) j 

= I dy + 

5 In the following, we denote the 4-momentum of this particle q = (qo,qr,q z ) and qr = |gr|. 

6 Defining x 1 = s^h/s, one indeed gets x-y = x-/(y, q y, Q, t ) where e = ±1 when the projectile (target) is 
the photon emitter. See the appendix A for details. 
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-3 -2 -1 0 1 2 3 4 J cms 

(a) proton on Pb 




-4 "-2 ' 0 ^ 2 ' 4 


(c) proton on H 

Figure 2: ^ for the proton run on Pb nucleon target (a), and for Pb run on H target (b), 
and (c) for p-p run. The flux is calculated for the specified Q 2 at qx = 0 (e.g. t = to). 
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where 4h^A)( y ) = s i h ( B , A ) from (A,B)(y)) and therefore 


^ = *£- ^( 4b ( 9)) + ^ d ^ hM) . (2.9, 

We therefore find it instructive to plot ^ (for fixed qr and Q 2 ) for different configura¬ 
tions: Fig. 2a shows the case of 7 TeV protons on lead, the lead being the photon emitter; 
Fig. 2b shows the case of 2.76 TeV lead on protons, the lead being the photon emitter * * 7 ; 
Fig. 2c shows the case of 7 TeV protons on proton (i.e. hydrogen), where both can be the 
photon emitter (note that the corresponding fluxes cannot simply be summed). 

In the first case (Fig. 2a), the flux is maximum in the forward region, which corresponds 
to a soft (coherent) emission by the lead target. In general, the particle tends to be 
emitted in the “hadron-receiver” region since the photon momentum is very small (here 
yreceiver=beam = 4.8). To be more precise, the boost - or rapidity difference Ay - between 
the photon-hadron and hadron-hadron cms, Ay'j^" = y^As 3 ~ Hems ' 3 , is simply related 
to the photon momentum fraction via 

= -4 lni ” (2-in) 

In order to produce a particle of mass Q at threshold 8 in a photon-hadron collision at y/s^h 
resulting from a hadron-hadron UPC at yJs^N-, one approximately has x 7 ~ Q 2 /snn- At 
AFTER@LHC, for Q 2 = 4 GeV 2 , x 1 = 3 x 10 -4 , which gives Ay^ ,h J l ’ 3 ~ —4. This explains 
the maximum at y cms ~ 4 of the solid line of Fig. 2a. 

If the emission is too soft, there is simply not enough energy to create a particle of 
a given mass, Q, - this explains why the curve for Q 2 = 4(40) GeV 2 drops sharply at 
y r ms — 4(3). In the second case (Fig. 2b), the flux is the highest in the opposite direction 
where soft photons are emitted by the lead projectile. Since the lead beam energy is lower 
{ybeam = 4.2) and the energy cut-off smaller, the particle is less backward than it is forward 
in the first case. In the third case (Fig. 2c), both proton can emit. The behaviour is similar 
to (a) and (b), but for a lower value due to the absence of the factor Z 2 . It is however 
harder in the forward (backward) region for a projectile (target) proton emitter - note 
that the flux in the tail is still nonzero at y down to ±3. This mirrors the possibility for 
harder photon (up to x 7 ~ 0.1) emission from a proton compared to a larger nucleus. The 
solid and dashed curves on Fig. 2c refer to two different fluxes: Drees-Zeppenfeld as in ep 
collisions (solid) and 6 m i n = 2 x R p (dashed). 

The advantage offered by the fixed-target mode is therefore obvious when the emitter 
is the projectile. In such a case, the large rapidity differences between the photon-hadron 
and hadron-hadron cms and that between the hadron-hadron cms and the laboratory frame 
nearly cancel. This results in the production of the particle at slightly positive rapidities 
which are easily covered by typical detectors. 

'In this case, as for (a), the photon emission by a lead ion can, in principle, be tagged with a neutron 

emission. In addition, the probability for the proton to be the emitter is, in practice, negligible ( Z 2 

suppressed at the same |i/ cm ,,|). We therefore do not show the curve for this possibility. 

8 i.e. at rest in the photon-hadron cms. 
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2.3.2 Double-photon case 


Obviously, one can also consider UPC where both colliding hadrons radiate a photon. This 
is expected to be the dominant reaction at work in dilepton production, AA' —> A£ + t~ A!, 
PbH—Tbf + r H, pPb—» p£ + £~ Pb or pp' —> p£ + £~p r , via 77 —£ + £~ (see Fig. 4a), that is 
the BH process. Combining the fluxes from both hadrons, one can derive a joint photon 
flux or 77 luminosity as function of their invariant mass, ^/s 77 or W. and rapidity Y. as it 
is usually done for gg luminosity at the LHC to discuss H° production rates, for instance. 
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10 1 
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0.001 


Figure 3: 


dL^y 


dWdY 


(and ) as a function of the rapidity, in the hadron-hadron cms 


W=m v 


(lower x-axis) and in the laboratory frame (upper x-axis) 


Fig. 3 shows the rapidity dependence of 77 luminosity in the hadron-hadron cms at a 
fixed W, obtained from of Eqs. 42, 43, 50, 51 and 52 of [27] with numerical integrations 
on the impact parameters. In the case of the BH process discussed in detail in the next 
section, the invariant mass of the pair is that one the dilepton. One observes that the 
maxima in the Pbp and pPb fluxes occur at respectively backward (forward) rapidities 
because the photon spectrum from the proton is harder than that from the ion. 

3 Lepton-pair production: energy, invariant-mass, rapidity and transverse- 
momentum dependencies 

We now discuss in more details the dominant process involved in lepton-pair production 
that is BH from 77 —»• £ + £~ (Fig. 4a) in hadron-hadron collisions. 

!, The reason why we took W = m rlc = 2.98 GeV for this example will become clear in the next section. 
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(a) BH (b) TCS with gluon (c) 7 - IP fusion (d) O - IP fusion 
GPDs 

Figure 4: Feynman graphs representing the scattering of two hadrons of momentum p\ 
and P 2 remaining intact with final state momenta p'i and p 2 producing a lepton pair via 
(a) the BH process by photon-photon fusion, (b) via TCS described in terms of a GPD 
(here the gluon one), (c) via photon-pomeron fusion represented by the minimal number 
of gluon exchanges and (d) via odderon-pomeron fusion also represented with the minimal 
number of gluon exchanges. 


System pp pYb Pb p PbPb 


y/swN [GeV] 

115 

115 

72 

72 


BW 

SL 

BW 

SL 

BW 

SL 

BW 

S 

dau(Q dQ GeV) [nb/GeV] 

0.14 

0.15 

200 

210 

77 

84 

7000 

7100 

t nb / Gev ] 

0.039 

0.038 

39 

45 

14 

16 

5500 

5600 

dau{Q= *Q 8 GeV) [nb/GeV] 

0.03 

0.031 

32 

34 

9.7 

11 

230 

250 

dauiQ=2 ZSi^ =0) i nb / GeV ] 

0.009 

0.009 

5.7 

6.5 

1.3 

1.6 

200 

210 


Table 2: BH differential cross section for fixed dilepton masses integrated (or not) on y^ms- 
“BW” denotes Breit-Wheeler and “SL” denotes Starlight. 

3.1 Total cross section and for fixed Q 

The total cross section for dimuon production by two heavy and charged particles is 
well known and can readily be computed from the analytical formula of Racah [28]. At 
y/s = 115 GeV, one gets 15.0 nb for pp collisions, to be compared with 16.5 nb with the 
Starlight Monte Carlo code 10 . The corresponding cross section for pPb is simply ob¬ 
tained by multiplying by Z 2 , i.e. 100 pb. In this case, it is assumed that the particle are 
point-like. Starlight rather gives 36 pb. 

We are however interested in the differential cross sections for particular values of the 
dilepton invariant mass, Q\ muon pairs produced at Q > c lm il are usually difficult to study 
since the muons have small transverse momenta. Such a cross section can easily be obtained 
by combining two EPA photon fluxes as done to obtain an ^ then by integrating over 

10 STARLIGHT website. 
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Figure 5: Cij x from protons and from Pb in pp, pPb, PbH and PbPb collisions (see 
Table 1) as a function of the invariant mass of the photon pair, W, at l^ ms =0 of the 77 
pair (or of the to-be produced dilepton). 


Y. One then convolve it with the Breit-Wheeler formula [29]: 


a 


.77 


(Q 2 ) = 


47 TOT 


8 mj 

~w 


1 6 mf\, Q + JQ 2 - 4m. 

1 ' In v 


Q 4 


2 m £ 


- \ 1 - 


4 mj 

~Q 2 


1 + 


Ynrq 

w 


(3.1) 


We have checked that we obtained the same result as Starlight for dimuon produc¬ 
tion, for instance for Q = (2, 2.98) GeV for pp and pPb collisions at y/s^N = 115 GeV, as 
well as for Pb p and PbPb collisions at y/sj\r n = 72 GeV, up to the uncertainty attached to 
the value taken for the nucleus radius, see the first and third rows of results on Table 2. 

One readily obtains the rapidity dependence of the differential cross section, by simply 
using jwdY I w=Q Y=y u ' a l rea dy discussed it on Fig. 3. The 77 luminosity at a fixed 
cms rapidity of the photon pair, y cms (or equivalently y^ms)) is plotted as a function of W 
on Fig. 5 namely at V- ms = 0 for pp, pPb, Pbp and PbPb collisions; it is multiplied by the 
corresponding hadron yearly luminosity denoted Cij. Just as for the rapidity integrated 
results, both methods agree for y^ s = 0 as it should be, see the second and fourth rows of 
results on Table 2. 

Knowing ^wdY 1 one can a ^ so obtain the production cross section for a scalar or tensor 
quarkonium, Q provided that we know its partial width into a photon pair, T 77 : 


da h n AhB 

~dy^~ 


7 7-^2 


’= 87^(2 J Q + 1) 


■ 77 


dL 


77 


2Mg dWdY 'W=Mq,y=Yq' 


(3.2) 


From r ^ 7 = 5.1 x 10 6 GeV and = 5.3 x 10 7 GeV [30], one gets a^ = 8.8 nb and 
= 3.2 nb. The rapidity dependence of this cross section is thus up to a constant factor 
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that of the joint flux at the corresponding 11 W (see Fig. 3). In particular, for pp (pPb) 
collisions at = 115 GeV, the r] c cross sections at y = 0 in the hadron-hadron cms 

are 0.5 pb (0.4 nb) and for Pbp (PbPb) collisions at snn = 72 GeV, 67 pb (11 nb). With 
the yearly luminosities in Table 1, one can respectively expect 10 , 1.8 x 10 6 , 74 and 80 r/ c 
per year. 

At this point, two remarks are in order. First, we stress that although, m Vc is above 
the energy “cut-off” showed in Table 1 for the systems pPb, Pb p and PbPb, the rates are 
nonzero using realistic photon fluxes. Second, our result for pp collisions is 2-3 times below 
that recently derived by Goncalves and Sauter [31] . We attribute this difference to the fact 
that they used the Drees and Zeppenfeld photon flux, which is precisely twice larger than 
that derived with b > 2 R p at x 7 ~ ^/m ric — 0.025 (see Fig. 1). As we noted above, the 
former should only strictly be used for ep collisions. Such a difference only arises at ’low’ 
hadron-hadron cms energies which correspond to x 7 on the order of 10 ” 2 and above. A 
measurement of the rj c or Xc 2 cross section in exclusive pp collisions at AFTER@LHC is 
therefore ideal to tell which choice is the most appropriate. 

3.2 Production at nonzero transverse momenta 

If we prefer to look at dileptons produced in reactions characterised by a momentum 
transfer squared of one emitter, |f|, up to 1 GeV 2 , which results in a measurable transverse 
momentum, qp. of the dilepton, it may be more suitable to treat the photon emission 
using a nucleon form factor 12 . By construction, the off-shellness of this photon cannot be 
neglected. This amounts to consider the process photoproduction, 7 p —>• £ + £~p, where the 
beam 13 photon flux is still treated in the EPA approximation and the effect of the form 
factor embedded in the photoproduction cross section. 



Figure 6: Definition of the kinematical variables in 7 (k)h(p) —>• £ + {l')£ (l)p(p') in partic¬ 
ular the angles 0 and attached to the £~. 

Following [ 8 ], this 7 h differential cross section for an unpolarised target to be convo¬ 
luted with the beam photon flux (provided that sin# m^/Q) reads 

11 Hence our choice of W = m^ c = 2.98 GeV in Fig. 5. 

12 In the nuclear case, such a configuration is admittedly much suppressed. 

lli The term “beam” may be improper in the fixed-target case since this photon can very well be emitted 
by a nucleon or nucleus in the target; this is particularly true for pPb collisions. 
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(3.3) 


aa ll _ ( U 

dQ 2 d.q^ d cos 0 dcj) cms 


a 3 

j em 


27rs 7 h 


11 + cos 2 0 
sin 2 0 


x 


(m - +(Fl(t)+ m)f 


where F\ 2 are respectively the Dirac and Pauli form factors evaluated at t = (p — p') 2 = 
(k — q ) 2 , the quantity e and procedure to obtain t(y^ ns ,q^,Q,e) and J (the Jacobian to 
change from t to q^) are explained in the appendix A. The 0 and (j> angles are defined 
as on Fig. 6. The apparent divergence at 0 —> 0 is regulated when the lepton mass is 
kept, just as in the Breit-Wheeler equation, Eq. (3.1), which is logarithmically divergent 
for mi —> 0. We note that 0 - the polar angle of the lepton - which is defined in the 
rest frame of the dilepton, can only be approximately related to the (pseudo)-rapidity of 
the lepton in the cms frame for ycms — 0 and qx -C Q. In such a case, the particular 
configuration (9 — 5 - 7t/ 2 corresponds to y^. ms — ?y^ ms ~ 0 which falls into the acceptance of 
a detector like LHCb and ALICE (with Pb beams). Using Eq. (2.9), one then obtains 
d&ie BU /(dQ 2 dq?p d cos 0 dcf> dy^ ms ). 


4 Timelike compton scattering or exclusive photoproduction of a dilep¬ 
ton 

The process of the lepton-pair production in ultraperipheral collisions may also be used to 
investigate hadron structure in terms of GPDs through the measurement of the contribution 
of the TCS process, Fig. 4b, to the cross section. Although the BH amplitude squared 
is much larger than the TCS one, it is possible to study the interference term between 
TCS and BH processes, which may be projected out through the analysis of the angular 
distribution of the produced leptons [8] and which depends on the GPDs. The interference 
term of the differential cross section is given by: 


da 


7 / 1 INT 

a 


oct, 


Vto - t y/1 - rj 2 


dQ 2 dt d cos 0 dcj) 4-7T s 2 h —tQ 

x Re 


COS (j) 


1 + cos 2 0 ' 
sin 0 


(4.1) 


Fi(t)H(ri,t) - v(Fi{t) + F 2 (t))U{r],t) - F 2 (t) £(rj,t) 


where we have neglected lepton mass and assumed that Q 2 ^$> t, M'jr. The variable p 
called skewedness is given by: 


Q 2 

2 s 7 h - Q 2 ’ 


(4.2) 


and to is maximal value of squared momentum transfer t reached at qx = 0, and is equal 
to = -4M 2 j)A up to corrections in 1 /Q 2 . The functions T-Ljp, t), t) and £(p, t) are the 
well known Compton form factors. These involve the GPD H, H and E (as defined in [15]) 
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respectively through a convolution with the hard-scattering kernels E computed at 
a given order in a s : 

{n,H,£}(ri,t)= (4.3) 

[ dx + T^ HSE} (x,r]){H 9 , H 9 , E 9 }{ X ,T],t) . 

The expression for the kernels T^ H E E -, at LO are given in the appendix B. The Next-to- 
Leading Order (NLO) hard-scattering kernel 11 , which we use here, can be found in [9, 32]. 

To advocate that the TCS measurement is feasible in the AFTER kinematics, we 
present phenomenological predictions making use of two GPD models through double dis¬ 
tributions [13]: the first is based on the Goloskokov-Kroll (GK) model based on fits of 
meson electroproduction data and the second is a model using MSTW8 parton distribu¬ 
tion function with a simple factorised t dependence, referred to as G-MSTW (for a detailed 
description see [10]). 



Figure 7: (left) <*^(66)/<^(66 = vr/2), (right) ^/^(69) for the model KG 
for Q 2 = 4GeV 2 , t = —O.lGeV 2 , and (j) = 0 for proton run on Pb, integrated over 
6 G (tt/4, 37t/4 ) and for y = 0 at NLO. 


Comparing the 4> dependence in Eq. (3.3) to that in Eq. (4.1), one directly sees that 
one can project out the interference term and eliminate the unwanted BH background by 
integrating the differential cross section over <f> with the weight cos cj). To avoid the limit 
6 0 where the BH signal gets very large, one can integrate over 8 symmetrically around 

7t/ 2 up to a value which depends on the size of the interference. Varying the limits of 
integration in 9 according to 


ir — 59 7T T 59\ 

2 ’ 2 ) ’ 


(4.4) 


results in the change of the strength of the interference BH-TCS with respect to the BH 
cross section shown on Fig. 7. On the left panel, we show the ratio of the the integral of 


14 As done in [9], we set = fx F = Q in the hard-scattering kernel and, as for other similar phenomeno¬ 
logical analyses, the GPDs are not evolved. In the case of the G-MSTW model, the GPDs can be made hf 
dependent through the evolution of the input PDFs. 
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the interference BH-TCS for a given 56 normalised to that for 59 = 7 r/ 2 . We see that for 
56 close to 7 r the ratio of the integrals of the interference stops increasing. In any case, 
in this limit, as we present on the right panel, the ratio of the integral of the interference 
term to the BH one gets very small because the BH cross section gets large 10 . In other 
words, there is no specific reason to take 56 much larger than 7 t/2 ; below this value the 
magnitude of the interference term is close to 10 % of the BH term. 

One can then combine the 7/1 cross section with the photon flux as a function of the 
dilepton rapidity at a fixed t to get the interference part of 5-fold hadron-hadron differential 
cross section. Defining the integration region as 


d b F = 



the BH cross section from a 7 TeV proton beam on a Pb target is 


a pPb = f (j 5 F _^BH_ 

BH Jn dQdq?pdyddd(f) 


1940 pb, 


(4.5) 


(4.6) 


which, for a luminosity of 0.16 fb 1 yr 1 , gives about 3 • 10 5 events per year. In the case 
of a H target, one has 


Q-PH = /' d 5p _^BH_ 

BH Jn dQdq^dydddcf) 


7.1 pb, 


(4.7) 


which, for a luminosity = 20 fb 1 yr 1 , gives 1.4 • 10 5 events per year for a 100 cm liquid- 
hydrogen target. Finally, for Pb on H, one has 


PbH _ 
^BH — 


= / dPF 


dciBH 


n dQdq^dydddcj) 


= 5500 pb, 


(4.8) 


which, for a luminosity = 11 nb _1 yr -1 , gives 6 • 10 3 events per year for a 100 cm liquid- 
hydrogen target. As aforementioned, we do not consider the case where the nucleus is 
emitter with a |t| of a few hundred MeV 2 which should be treated with nucleus form factors 
and nuclear GPDs. With a magnitude of 10 % for the interference term, the azimuthal 
modulation should be observable in the 3 cases. 

On the Fig. 8 , we show the magnitude of the 3 terms as function of the rapidity as 
well as the relative magnitude of the interference term with respect to the BH one for three 
cases of collisions (a) proton beam on a Pb target, (b) Pb beam on a H target, (c) proton 
beam on a H target. On the left panels, we present the rapidity dependence of the BH, 
TCS and BH-TCS terms of differential cross sections evaluated at Q 2 = 4GeV 2 , 

t = — 0.35GeV 2 and <f> = 0 and integrated over 6 G (7t/4, 37t/4) for the GK model. On 
the right panels, we show the ratio of the interference signal to the BH for the GK and 
G-MSTW model. We see that for all cases the signal to the background ratio is around 
10-15%. 

15 In one neglects the lepton, it even diverges. 
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(a) 7 TeV p on Pb target 



(b) 2.76 TeV Pb on H target 




(c) 7 TeV p on H target 

Figure 8: (Left) Differential cross section dydt c ^ d Q 2 for the KG model for Q 2 = 4GeV 2 , 
t = —0.35 GeV 2 and <j> = 0 integrated over 6 6 (vr/4, 37 t/4 ). Dotted line : B-H, dashed line 
: interference term, solid line: TCS. (Right) Ratio of the interference to BH differential 
cross section dydt d d ^ d Q 2 calculated for Q 2 = 4 GeV 2 , t = —0.35 GeV 2 and </> = 0 for the GK 
(dashed) and G-MSTW (solid) models at NLO.From (a) to (c) pPb, PbH and pH cases. 
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5 Exclusive lepton-pair hadroproduction via odderon-pomeron fusion 

Throughout this work, we have based our discussion on the idea that the (theoretical and 
experimental) requirements for selecting UPC collisions was de facto preventing hadronic 
interactions to take place or at least that photon-induced processes would be dominant. 
This is a rather safe assumption in nucleus-nucleus collisions where the coherent photon 
fluxes are enhanced by a factor Z 2 for each nucleus, whereas double-pomeron-induced 
reactions only scale like A 1 / 3 . It also seems a reasonable assumption for pA collisions that 
electromagnetic interactions dominate. 

In the case of pp collisions, the Z 2 factor is absent and the UPC requirements are also 
very similar, if not identical, in practice, to that for an exclusive or diffractive process. For 
instance, in order to impose b > 1 /R p , one can impose that the proton expected to emit 
the photon is only deflected by a Apt of say maximum 100 MeV, which corresponds to 
b > 2 fm. Such a requirement however may not be sufficient to prevent any scattering by 
the exchange of a pomeron or an odderon. In fact, the photon-pomeron process depicted 
on Fig. 4c should in principle be encompassed in the GPD description at high energies. 

The production of a lepton pair in UPC which we discussed so far should thus be 
confronted to a potentially competing process leading to the same final state, in which the 
virtual photon decaying into lepton pair is produced in an Odderon-Pomeron (OP) fusion 
as depicted in Fig. 4d. This is similar to the situation occurring in the case of the central 
exclusive hadroproduction of a J/if studied in [33]. In this case, the J/if can be produced 
either by the photon-Pomeron fusion - that is the one expected for UPC - or by the OP 
fusion. In fact, the cross sections for both processes are similar at the Tevatron. Owing 
to the increase of the gap survival probability for lower energies, it is reasonable to believe 
that, at AFTER@LHC energies, i.e. around y/sftN = 115 GeV, exclusive J/if production 
in pp collisions would be dominated by OP fusion. 

Both production mechanisms are however characterized by a different J/if transverse- 
momentum dependence [33]. In particular, if one imposes that the final state protons have 
|t| > 0.25 GeV 2 , the OP contribution can be enhanced ten times with respect to the photon 
induced one. 

In the case of lepton-pair production in the conditions of UPC, a similar situation will 
likely occur. The experimental study of such a process with specific cuts on t may thus 
provide an unexpected path to the Odderon discovery via interference with the pure QED 
BH process in the same way as the TCS signal would be extracted in the region where it 
dominates. Despite its interest, the evaluation of the corresponding cross section along the 
line of [33] is far beyond the scope of the present work. 

Similarly, the odderon could also contribute to the rj c production as recently discussed 
in [31]. rj c can be produced by photon-odderon fusion due its different quantum number 
with respect to the J/if- In the pp case, the cross section we obtained from 77 fusion is 
slightly larger that from photon-odderon fusion obtained in [31], both on the order of a 
picobarn or less. In the pPb and Pbp cases, the photon can be radiated by the ion and 
the pomeron by the proton. As we noted, although m ^ is in principle above the so-called 
energy “cut-off”, the photon fluxes, even accounting from the minimum impact parameter 
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for a proton-lead UPC, is not zero. In the Pbp case, the photon-odderon induced cross 
section ranges from 30 up to 360 pb whereas the 7 7 induced one is on the order of 440 pb. 
As for the dilepton case, a study of the transverse-momentum dependence should be able 
to discriminate between both processes. 

6 Conclusion 

We have theoretically investigated the feasibility of accessing the lepton-pair production 
in ultraperipheral collisions at the proposed fixed-target experiment AFTER@LHC, which 
takes advantages of the multi-TeV proton and ion beams of the LHC. To this aim, we 
have first estimated the magnitude of the cross section for lepton-pair production from 
the fusion of two quasi-real photons emitted by the quasi-grazing hadrons. This purely 
electromagnetic BH process can serve as an important tool for the determination of the 
luminosities with nucleon or ion beams but it can also be used for an experimental veri¬ 
fication of the validity of the effective-photon approximation usually applied to estimate 
the flux of quasi-real photons emitted by these relativistically moving charges. Lepton-pair 
production also gives access to the proton GPDs via the TCS process. Another way to 
probe the photon flux is to measure r) c production for which the production rate in pp 
collisions at AFTER@LHC does depend on the minimal impact parameter used for the 
UPC. 

The predictions which we obtained for the cross section for BH - using specific cuts 
relevant for the GPD extraction - are on the order of few thousand of pb for the pPb and 
Pb p cases and a slightly less than 10 pb for the pp case and we confirm the dominance of 
BH over TCS. This dominance can partially be overcome by studying of the interference 
- also sensitive on the GPDs - between TCS and BH which we evaluated at NLO. With 
specific cuts on the lepton polar angle, the ratio of this interference over the BH amplitude 
squared is on order of 10% with two models of GPDs, i.e. GK and G-MSTW. These 
are quite promising values giving hope for the extraction of this interference by means 
of the analysis of the azimuthal distribution of the produced leptons. Studying TCS in 
ultraperipheral collisions at a fixed-target experiment can also give us opportunity to study 
target polarization asymmetries, which are an useful tool to extract further information on 
GPDs [34], 

We have also derived cross sections for r/ c production by photon-pair fusion, which 
happens, in particular in this energy range, to be sensitive on the method used to compute 
the photon flux in the pp case. We found out that 10 4 r/ c can be produced per year in 
UPCs with AFTER@LHC. 

Finally, we discussed possible competing hadronic processes via pomeron or odderon 
exchanges which could interestingly be separated out by a careful analysis of the transverse- 
momentum dependence of the produced particles. 

In conclusion, AFTER@LHC offers a realistic possibility to study lepton-pair produc¬ 
tion in ultraperipheral collisions which opens the path to investigate features of the partonic 
structure of hadrons which are complementary to those studied with lepton beams. The 
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use of hadron beams may, for instance, offer the opportunity study to odderon-sensitive 
reactions. 
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A Kinematics 

We denote the momenta of incoming nucleons (in the nucleon-nucleon cms) as: 

PA = (1,0,0, a), 

PB = (!, 0, 0, —a), (A.l) 

where the A is a nucleon from the beam, B is a nucleon from target and a = \J 1 — 

The Weizsacker - Williams photon is emitted from a beam (e = —1) or a target nucleon 
(e = +1), and its momentum is given by: 

k = x^( 1,0,0,-e). (A.2) 
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Momentum of outgoing lepton pair (or outgoing virtual photon decaying into heavy lepton 
pair) reads: 

q = ( qo,qT,q z ) = (m T cosh y u ,q T ,m T sink y u ), (A.3) 

where mx = \J qr 2 + Q 2 and y u is the lepton-pair rapidity which can be expressed as 

sx 2 (a + 1) 

where: 

t = (k — q) 2 = Q 2 — mT'/sx 1 e eyU . (A. 5) 

Inverting Eq. (A. 4), we can express as a function of y . 


(A.4) 


I e log ( _ (Q 2 ~t)(a + 1) 

2 yQ 2 (a — 1) — t(a — 1 — 2x 7 ) 


x^{y U ,t,Q,e) = (A.6) 

—2 1 + ^4t 2 — 4s(Q 2 — i)(o! + l)[(a — 1) — (a + l)e -2e ^] 

2s(a +1) ' 

or combining Eq. (A. 5) with Eq. (A. 7) we easily get: 


t = t{y ee ,q^,Q,e), 

X", = x^{y u ,qT,Q,e), 

_ dt 

dq^ 


(A.7) 

(A.8) 

(A.9) 


B Compton form factors and generalised parton distributions 


In this appendix, we give the expressions of the LO hard-scattering kernel T % appearing in 
the expression of the Compton form factor Ti, 7-L and £. At Born order, the hard-scattering 
kernel associated to the quark GPDs are given by 


T q H (x,y) 
T]j(x, rj) 

T q E (x,v) 



r] — x — ie —y + x — xe , 


(B.l) 


and those associated to the gluon GPDs are zero. The NLO hard-scattering kernels can 
be found in [9, 32], 


- 21 - 










